Prom particle in a box to PT-symmetric systems via isospectral deformation 
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A family of PT-symmetric complex potentials are obtained which is isospectral to free particle in 
an infinite complex box in one dimension (1-D). These are generalizations to the cosec 2 (x) potential, 
isospectral to particle in a real infinite box. In the complex plane, the infinite box is extended 
parallel to the real axis having a real width, which is found to be an integral multiple of a constant 
quantum factor, arising due to boundary conditions necessary for maintaining the PT-symmetry 
of the superpartner. As the spectra of the particle in a box is still real, it necessarily picks out the 
unbroken PT-sector of its superpartner, thereby invoking a close relation between PT-symmetry 
and SUSY for this case. 
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Introduction 
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Supersymmetric quantum mechanics [l|, Q provides an elegant connection between two different systems with 
Hamiltonians in factorized form 0] . The additional property of shape- invariance Q ensures solvability of the potential 
by relating it to a family of systems with successive isospectrality Interesting phenomena, like refiectionlessness 
potentials Q can be explained to be due to isospectrality of the system Q with free particle. On the other hand, 
PT-symmetric systems @ have been under consideration owing to the realness of the spectra belonging to non- 
Hermitian complex Hamiltonians (potentials). These turn into complex-conjugate pairs of energy for a different 
range of parameters interpreted as spontaneous breaking of PT-symmetry. For unbroken PT-symmetry, the realness 
of the spectra necessarily demands generalization of the usual Dirac-von Neumann scalar product [7[, generically 
f*"*. ■ unique to the particular system, in the same spirit of the pseudo-Hcrmitian systems. It should be emphasized that 
the equivalence of PT-symmetry and pseudo-Hermiticity is yet to be established § . Certain class of PT-symmetric 
potentials display isospectrality to Hermitian systems . This is expected to hold in the same parameter range, where 
t-h ' PT-symmetry is unbroken. The necessary realness of both spectra is ensured by supersymmetry. Thus, a co-existence 
t-H . of both the symmetries is observed. Recently, supersymmetric structure of certain PT-symmetric systems has been 
considered to analy_ze the structure of the corresponding Hilbert spaces, along with the parametric correspondence 
between the two [lOj . Subsequently, scattering properties of PT-symmetric systems were qualitatively studied to 
realize unique boundary conditions arising for such systems [Tlj j . in accord to the experimental results 

In this paper, we will analyze and generalize the simplest quantum mechanical system, the one dimensional infinite 
well, to generate a whole family of PT-symmetric potentials through isospectral deformation 0. Strict boundary 
conditions necessarily require complexification of the 1-D space, thereby resulting in multiple choices for the parity 
transformation. We stick to the definition of parity that alters the sign of only the real part of the co-ordinate, in 
order to stay close to physical realization. It follows that the width of the well, which has to be real for the realness 
of the physical spectrum, cannot be arbitrary. Instead, it has to be an integral multiple of a constant, manifestly 
quantum in nature. PT-symmetric systems, under constant imaginary shift of coordinate, have been studied |13j for 



their supposed '?7-pseudo-Hcrmiticity' 14} . Hhere r\ is the proposed 'norm operator' for the pseudo-Hcrmitian systems 
[ToT ]. In the present case, we utilize 'general isospectral deformation' to arrive at complex shift in coordinate of the 
'real' system, obtaining a map to a PT-symmetric system, defined in real space with imaginary coordinate acting as 
a continuous parameter, subjected to the physical boundary conditions. A scheme of extending the isospectral family 
of potentials to include non-Hermitian systems has been realized. 
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The paper has been organized as follows. Section I contains a brief summary of isospectral deformation leading to 
the solutions of the Bernoulli's equation. Section II deals with the family of equivalent superpotentials corresponding to 
1-D particle in an infinite box and that the resulting isospectral family of potentials necessarily require complexification 
of space and discretization of the well width. Section III contains a generalized treatment of 1-D infinite well in complex 
plane, emphasizing on the importance of boundary conditions necessary for physical results. Section IV analyzes the 
■pT-symmctric potentials, isospectrally constructed from the generalized particle in a box, and also the parametric 
conditions leading to spontaneous breaking of "PT-symmetry, thus shading light on its effect on the isospectrality. In 
conclusion, the accumulated understanding of these generalized systems will be discussed, along with the prospects 
for future works. 



I. ISOSPECTRAL DEFORMATION 



Given a potential leading to a factorizable Hamiltonian, the corresponding superpotential is not unique. However, 
different superpotentials correspond to different isospectral partners of the unchanged potential in general. Such 
modification of the superpotential is known as the isospectral deformation, and evidently in this way, a family of 
superpartners can be constructed for a given potential. An isospectral pair of potential is conventionally written as, 

V-(x) = W 2 (x) - — =W'{x) 
\j2m 

V + (x) = W 2 (x) + -jLW(x), (1) 

which have exactly the same energy eigenvalues, except for the ground energy of V-(x) with no mapping into V + (x) 
[i.e. = E~ +1 ], but different cigenfunctions which are algebraically related as [l|, 

il>+(x) = [E- +1 ]-iKip- +1 or VWi(z) - [E+^K^-. 

Here, W(x) is the superpotential and W'(x) is its first derivative with respect to x. We invoke the isospectral 
deformation W(x) — > W(x) + g(x), g(x) being an arbitrary function restricted only by the boundary condition 
of the unchanged system, chosen for convenience to be that for V-{x), and normalizability of the corresponding 
wavefunctions for the new V+(x). Then, from EqQ] for unchanged V-(x) we obtain: 



\x) + 2W{x)g(x) 



1 2m 



g'(x) = 0, 



(2) 



with a new superpartner, 



V + (x) = (W(x)+g(x)f 



^2mdx 



(W(x)+g(x)). 



(3) 



The physically acceptable solution g(x) of the Bernoulli's equation (Eqj2j), for a given W(x), is determined by the 
normalizability of the eigenfunctions corresponding to V+(x). 



II. ISOSPECTRALITY OF THE FREE PARTICLE IN INFINITE BOX 



The free particle corresponds to a constant superpotential W(x) = A. Presence of local boundary conditions 
generally modifies it to a local function, as in the case of 1-D infinite box, reflected in the expression of the ground 
state in terms of the superpotential, ipo(x) = exp [— J W(x')dx'] . The constant potential A 2 can always be absorbed 
in defining the origin of energy. The corresponding Bernoulli's equation reads as, 



g 2 (x)+2Ag(x) 



'2m 



g'(x) - o, 



with the solution, 
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g(x) = -2A- 



,2A(- 



fiA(- 



■-+c) 



(4) 



with c being the constant of integration. The fact that V~(x) is a constant, leads to the Riccati equation: W 2 (x) — 
W (x) = constant for the superpotential, having unique solution modulo an additive constant. Though it is 



/2 m 



trivial in this particular case, it ensures uniqueness of 9 {x) too. The new superpotential -Acoth (a^x + Ac 
which will hereon be called W(x), generically marks the presence of local boundaries and leads to a new superpartner, 



V+{x) = 2A 2 sinir 2 A- 



Ac 



A 2 



(5) 



It is to be emphasized that the functional form of the potential (V-(x)) by itself does not contain information of 
the boundary conditions, but the superpotential does, as it determines the ground sate eigenfunction. Thus, through 
isospectral deformation, only the value of V_(x) is unchanged, not the system boundaries. 
The ground state wave function for V-{x) in terms of the new superpotential is, 



^{x)=N[-e 



l -x+2Ac 



(6) 



which must be same as that for W(x) = A, i.e., sinusoidal with nodes only at the boundaries. This necessarily requires 
that A = ia, where aeR, yielding, 



V+(x) = 2a 2 esc 2 



(7) 



The knowledge of the ground state of V-(x) along with EqJS] ensures the trivial physical choice of c = 0. This is 
the previously mentioned boundary value restriction from a potential (V-(x)) to its superpartner (V + (x)), as evident 
from Eq[7l However, we are motivated here to construct the most general form of V+(x), thus to generalize V-(x) if 
necessary. Keeping this in mind, we consider the integration constant c to be complex in general. Then, subjected to 
the boundary conditions that the eigenfunction vanishes at X12, one obtains from EqJHl 



J2iac 



-2ia± 



-Xl,2 



(8) 



thereby necessitating x\.2 to be complex. The above equation also ensures the difference between x\ and xi to be 
j J^ m — , with j being an integer. Therefore, we end up generalizing the 1-D infinite potential well with arbitrary width 
to be an infinite 1-D potential well in the complex x plane, placed parallel to the real axes, with real but integral 
'quantum' width. That this choice is the most general one, is explicated in the appendix. Thus, V— (x) can correspond 
to different infinite wells subjected to different combinations of j and a. The eigenvalues, however, are still real as 
they depend only on x\ — X2, which is real, but can have different values now. Hence, we ended up isospectrally 
relating a family of 1-D infinite wells to a family of complex potentials, which will be discussed next. 
For mathematical and intuitive convenience, we redefine the variable(s) as, 

p = a — - — x re + ac r 



-x m + ac m , with 



(9) 



Then, one can re- write x — ° \(p — ac re ) + i{q — 



in terms of the new variables. Now, the potential in Eq(7] 



can be written as V+(x) = 2a 2 cosec 2 {p + iq) — a 2 which can be re-organized to become, 



V+(x) 



^ sin p cosh q — cos 2 p sinh q 
(cos(2p) - cosh(2g)) 2 



sin(2p) sinh(2q) _^ 
(cos(2p) — cosh(2<7)) 2 



(10) 
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(a) (b) 



FIG. 1: Plots of (a) Me(V+(x)) and (b) Im(V+(x)) vs. x re modulo an additive constant c re (or p). 

with real and imaginary parts well-separated. Fig [T] shows the plots of V+ e (x) and V] r m (x) as a function of x re , 
showing the former to be even and the later to be odd in x re . This ensures PT-symmetry of V+(x), provided we 
define parity transformation in the complex x- plane to be x re — > —x re ,Xi m — ¥ Xi m . This is sensible also from the 
point of view that the generalized potential wells are actually manifested through the shift of the origin, as far as the 
invoked complexity of the coordinate x is concerned. The variable Xi m acts as a continuous parameter extending from 
+00 to —00, same as a and c, while the variable x re is restricted by boundary conditions for the well. 

Hamiltonians in complex coordinate plane have been extensively studied in classical domain (l6j , wherein multiple 
trajectories, along with the shifted 'usual' real space trajectories (straight lines), correspond to the same probability 
in corresponding quantum systems [13] • These have different sets of 'classical turning points' with the possibility of 
'classical tunneling Here, we stick to the shifted 'usual' trajectory for the infinite well, which is fixed parallel to 

the x re axis by the physical boundary conditions to be explicated in the next section. Subsequently, V+(x) is confined 
in real space only, with Xi m serving as a continuous parameter. 



III. THE GENERAL 1-D INFINITE BOX IN COMPLEX SPACE 



In order to obtain an isospectral map to 'PT-symmetric system of V + {x) in real space, by generalizing a Hamiltonian 
system (the 1-D infinite box) to the complex plane, we need Xi m to be at best a continuous parameter, whereas x re is 
the direction along which parity transformation has to be executed. Here, one expects the motion of the particle in 
the infinite well to remain confined to the real axis, for the spectra of the system to be real, requiring the quantized 
widths of the well to be real. This requirement makes the treatment fundamentally different from complex canonical 
transformation in the phase space, utilized to study PT-symmetric systems earlier [1 9j j . obtaining a 'new' Hamiltonian 
with real spectra. Here, we do a position-coordinate or 'point' transformation by taking the Hcrmitian system to 
complex plane, and then obtain PT-symmetric systems from there through isospectral deformation. 

We begin by considering the 1-D infinite box to be oriented in the complex plane [x re , Xi m ] in arbitrary direction 
instead of parallel to the x re -axis. Consequently, the generalized boundary conditions for the eigenfunction will be 
ipnix^e , x \m) = = ^nO^re , x \m)i with coordinates in the bracket representing the infinite boundary. Now, the 2-D 
Laplacian, V 2 = -Jp 8 d . A is defined in the complex coordinate plane, yielding the Schrodinger eigenvalue equation, 



q2 q2 \ *2rnE 

ip(x re ,x im ) = -K 2 ip(x re ,x im ), where K 2 = — . (11) 



ox re ox im J H 
On separation of variables as %jj{x rei Xi m ) — lZ(x re )I(xi m ), one obtains 



—1Z" + K 2 = K 2 = —I", (12) 
K 1 

where K 2 is a constant and the prime denotes derivative with respect to the respective variables. The individual 
solutions are: 
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lZ(x re ) = Aex.p(ikx re ) + B exp(— \kx re ), and X(xj m ) — C exp(kxi m ) + D exp(— kxi m ), (13) 

Where A, B, C and D are complex numbers and K 2 = K 2 — K 2 . We identify that K = K re and K = Ki m , the 
momenta along respective axes. 

The second one of EqfT^lis what one obtains considering a 1-D free particle in imaginary space. In that case, K 
is real as the space is considered to be imaginary. If one considers the variable of differentiation to be ixi m , instead 
of x, m , we end up with imaginary momentum. These are two equivalent pictures of a 1-D free particle in imaginary 
space, which can also be seen as a 1-D particle moving in a constant potential V > E, E being the total energy of 
the particle, known to have exponentially decaying solutions. 

On plugging in the boundary conditions ipn {%re , i|„ ) = = ip n (i™ , xf^ ) in Eq. [T5J we obtain, 

Aexp(ikx^) + B exp(-ikx^) = = Aexp{ikx {2 }) + B exp(-ikx<$) 
and C exp(kx[ll) + D exp(-kx^) = = C exp(kx^) + D exp(-kx^l) , (14) 

yielding, 

exp |lK"(x^ — x£^)j = 1 — exp(i2n7r) and exp |-K"(x-j^ — x-^) j = 1 = cxp(i2m7r), n,mel. (15) 

As K, K, x re and Xi m are all real, the only sensible conclusion is m — 0, implying x£^ — x\^. Thus the box is parallel 
to the real axis and the particle is constrained to move in 1-D. The identification K — Ki mi owing to the explicit 
consideration x, m eR, ensures the 'physicality' of the Hermitian system (1-D infinite well) that it has real spectra with 
the wave- function decaying exponentially along Xj m , with proper choice of coefficients C and D in Eq ll3l This is a 
general result, owing to both the imaginary nature of one variable and the infinite boundary. 

The above problem can equivalently be cast as one in the 5*0(1, 1) space (x, y) (the 1 + 1 Minkowski space), where 
the conserved scalar is x 2 — y 2 . The corresponding conserved scalar in the Fourier space is K 2 — K 2 , with the metric 
g^v = Diag(l, — 1) or = Diag(l, i), where the space is essentially considered to be real. We adopted the notation 
of the equivalent metric = (1, 1) (50(2)), with explicit consideration that the space is complex. This leads to the 
conserved scalar x 2 — y 2 same as that in the 1 + 1 Minkowski space, but now that in the Fourier /momentum space is 
K 2 + K 2 . Thus under the 50(2) metric, either space (1-D component) can be real or the corresponding momenta, 
but not both, if the complete space is complex. The essence to carry forward here is that in the present problem, 
the complexness of space is not avoidable as the 50(2) metric is unique, and hence the imaginary component of the 
space is not an artifact, removable by choice of metric, unlike the 50(1, 1) case. 



IV. THE GENERALIZED FAMILY OF PT-SYMMETRIC POTENTIALS 

The spectrum of V-(x) is still real, which can be verified in the usual manner by considering the general plane 
wave solution Ae lkx + Be~~ lkx to vanish at complex values 2+2, leading to k(x2 — x\) = mr with (x2 — x%) — L being 
real. This, through isospectrality, necessarily demands the spectrum of V+{x) to be real. But being 'PT-symmetric, 
V+{x) also has complex-conjugate spectra when the 'PT'-symmetry breaks spontaneously. Thus the isospectrality with 
the complex infinite well should persist only till the 'PT-symmetry is preserved. To analyze the complete spectra of 
V+ (x) , we seek to solve the eigenvalue problem for H + — — j— -4-z + V+ (x) . For convenience, we consider the form in 
Eql3 leading to the eigenfunctions [l|: 

+ / I \j2m I I \J1vn I / \J1vn \ 

iPn\ x ) ^ \ n + 1) cos { ( n + 2) (a — - — x + ac) > — sin < (n + l)(a — - — x + ac) > esc I a — - — x + ac I , (16) 

modulo a normalization constant. On considering Eq. [9l the above equation can be rewritten as: 

ipn( x ) ^ ( n + 1) cos{(n + 2){jp + iq)} — sin{(n+ l)(p + ig)}csc(p+ iq) 

oc (n + 1) [cos{(n + 2)p} cosh{(n + 2)q} — i sin{(n + 2)p} sinh{(n + 2)q}] 

— [sin{(n + l)p} cosh{(n + l)q} + i cos{(n + l)p} sinh{(n + l)q}] [sm(p) cosh(q) — i cos(p) sinh(g)] 1 (17) 
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(a) (b) 
FIG. 2: Plots of (a) M.e(V+(x)) and (b) Im(V+(x)) vs. Xi m modulo an additive constant Ci m (or q). 

Clearly, these eigenfunctions are PT-symmetric with respect to x re modulo a shift ac re of the origin. This is in accord 
with the property of isospectral PT-symmetric potentials that unique superpotential leads to unbroken PT-symmetric 
sector, as observed in [Toj . unless any further symmetry is involved [20]. 

From Eqs [POl and [171 it is also evident that both V+(x) and tp„(x) remain PT-symmetric subjected to the 
alternative choice of parity transformation: x re — > x re ,Xi m — > —Xi m , with x re acting as a continuous parameter 
(Fig [5]). However, this will not correspond to a particle in an infinite box directly, as such a system cannot exist 
parallel to the imaginary axis. Therefore, no alteration of the parameters can lead to a spontaneously broken PT- 
symmetric regime at this stage. Naturally, there exists a more general form of V+( x), which reduces to that in EqllOl 
under certain parametric condition, preserving PT-symmetry spontaneously [lOj . To obtain the same, instead of 

—A coth ( A ^™ x + Ac) , we propose a superpotential: 



W{x) 



-a cot a- 



'2m 



ac 



LB. 



-L4 and aeM, BeZ, 



(18) 



leading to, 



V + {x) = {a 2 + aa) esc 2 \^a ^ x + acj - i2aBcot + ac J _ ^ + B ^ ' '- 19 -' 

The potential in EqjT] corresponds to a particular choice a = a and B = 0. For a ^ a and B = 0, the system will 
not map to a free particle (in or out of the box!) as V-(x) — (a 2 — aa) esc 2 (a ^f™ x + acj — a 2 with ground state 

wave function ipQ (x) oc |sin (^a ^™ x + acj j . It is evident that V±(x) are shape-invariant P, [J |3| and hence, 

can be solved exactly using the SUSY algebra. The spectrum turns out to be E n = a 2 — (a — na) 2 p], [Hi, which is real. 
Therefore, the TT-symmetry is still preserved. On the other hand, for a ^ a and B ^ 0, it is known [10] that there 
exists a spontaneously broken PT-symmetric phase, with complex-conjugate pairs of discrete, yet finite eigenvalues, 
even for xel&. 

Therefore, in the complex space, the free particle in infinite 1-D box isospectrally maps to a restricted portion 
(a = a, B = 0) of the parametric domain for spontaneously preserved PT-symmetry of the potential in Eq ll9l Such 
results were known earlier [2l| for specific PT-symmetric potentials obeying SUSY-QM. Isospectral deformation plays 
a distinct role in the present case, by bringing out the most general form of the superpotential corresponding to a 
constant potential. Generic complexification yielded the PT-symmetric potential, isospectral to the free particle in a 
real 1-D infinite box, embedded in the complex space. 



Conclusions 



The fact that the real spectrum of a PT-symmetric potential can be mapped to that of a much simpler potential 
through suitable generalization of the space enables us to intuitively picture why the spectrum of such a potential is 
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real, even without Hermiticity, due to isospectrality to a Hcrmitian system. Here, we have extended the domain of 
such analogy through generic complcxification of the space. This might be helpful to visualize the real spectrum of 
other 'PT-symmetric potentials. This could serve as an alternative to constructing a new scalar product to explain the 
real spectrum of these complex potentials, by attributing the realness of the spectrum to parametric isospectrality. 
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